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Abstract 

We provide accurate Monte Carlo results for the free energy of in- 
terfaces with periodic boundary conditions in the 3D Ising model. We 
study a large range of inverse temperatures, allowing to control correc- 
tions to scaling. In addition to square interfaces, we study rectangular 
interfaces for a large range of aspect ratios u — L1/L2. Our numeri- 
cal results are compared with predictions of effective interface models. 
This comparison verifies clearly the effective Nambu-Goto model up 
to two-loop order. Our data also allow us to obtain the estimates 
r,/v^ = 1.235(2), TOo++/v^ = 3.037(16) and R+ = /^cto = 0.387(2), 
which are more precise than previous ones. 



1 Introduction 



Interfaces play an important role in various fields of natural sciences. In soft 
condensed matter physics, in chemistry and in biology, interfaces separating 
two different media, for instance two different magnetization domains, or 
two different fluids, or a fluid and its vapour, are studied. The properties of 
such interfaces might be described by a unique effective model such as the 
capillary wave model [1]. 

Our motivation to study interfaces originates from the theory of high 
energy physics. An interface with given boundary conditions can be asso- 
ciated with the world-sheet of a fluctuating flux tube in the confinement 
regime of a gauge theory. For intermediate and long distances between the 
sources, the relevant degrees of freedom for a system of confined quarks are 
supposed to be independent of the short distance gauge interaction, and 
might be modelled by string fluctuations {effective string picture). 

The simplest set-up for a numerical study of interfaces is provided by 
the Ising spin model on a simple cubic lattice. Its duality with respect to 
the Z2 gauge model [2] maps the ordered phase to the confined regime. 

The classical Hamiltonian of the Ising spin model reads: 

H{{J},{h},{s}) = -^Jf^^y^SxSy - ^ha^Sx , Sj,g{1,-1}, (1) 

(xy) X 

where x = {xo,xi,X2) is a site of the lattice, and (xy) denotes a pair of 
nearest neighbours on the lattice. Here and in the following, the lattice 
spacing a is set to 1, and we shall always consider the case of a vanishing 
external field hx = 0, Vx. The site coordinates run over < Xj < Lj — 1, 
where i e {0, 1, 2} label the three directions. 

In the case of periodic boundary conditions we take J{xy) = 1 for all 
links {xy). Anti-periodic boundary conditions, say, in the direction 0, can be 
implemented imposing J(xy) = — 1 if x = (Lq — 1, xi, X2) and j/ = (0, xi, X2), 
and J{xy) = 1 otherwise. 

The partition function is obtained as the sum over all configurations {s} 
of the Boltzmann factor: 

Z{jy{P) = J2^w{-m{J},{s})) , (2) 
{4 

where jS = 1/ (A'bT) is the inverse of the temperature of the three-dimensional 
classical spin model. 
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The goal of our work is to study an interface between the phases of 
positive and negative magnetization in the low-temperature regime of the 
spin model — which corresponds to the confining regime of the gauge theory. 

Such an interface can be forced into the system by appropriate boundary 
conditions. For instance, one could constrain the spins at xq = to take the 
value —1 and those at xq = Lq — I to +1; here, however, we use anti-periodic 
boundary conditions in 0-direction, because the finite Lq effects are smaller 
and better understood than for Dirichlet boundary conditions. 

In recent works [3-12] we studied interfaces with Dirichlet boundary 
conditions in one direction and periodic boundary conditions in the other 
direction: via duality, this corresponds to a Polyakov loop correlator in the 
gauge model. 

The comparison with the Nambu-Goto effective string model resulted in 
unexpected discrepancies at subleading orders. While finite L2 corrections, 
in the direction with periodic boundary conditions, are described well by the 
effective theory, the finite Li corrections, in the direction of the Dirichlet 
boundary conditions, show unexpected deviations. 

In order to further investigate this issue, we pick up again the work on 
interfaces with periodic boundary conditions in both directions. In [13- 
15] such a comparison had been performed for square interfaces Li = L2; 
in these studies, the numerical values of the interface tension were taken 
from [7]. 

In the present work, our results for the interface free energy allow for an 
independent determination of the interface tension, which is computed in 
technically quite a different way with respect to [7]; the consistency of the 
two results provides a non-trivial check of their validity. 

We obtain results for a large range of the inverse temperature f3, allowing 
to study possible scaling corrections. Furthermore, we also compute the 
interface free energy for Li / L2 for a large range ofu = L1/L2: this enables 
us to compare with the non-trivial dependence on n, which is predicted by 
the effective interface models. 

Finally, the results for the interface tension are also used in combination 
with a series analysis of the second moment correlation length in the high 
temperature phase. This yields a precise estimate of the universal amplitude 
ratio: 

^+ = /2nd, + <70 , (3) 

where the amplitudes are defined by cr ~ ao{—t)^^ and ^2nd — /2nd, Here, 
a is the interface tension, ^2nd the second moment correlation length in the 
high temperature phase, t = {T — T^)/T^ the reduced temperature and i^, 
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II = 2v the critical exponents of the correlation length and the interface 
tension, respectively. The result for can be compared e.g. with results 
obtained from experiments on binary mixtures. 
We also update the estimate for: 



where now the error is dominated by the estimate of the mass ?tio++ of 
the 0^"*" glueball. Note that under duality the interface tension of the Ising 
spin model is equal to the string tension of the Z2 gauge model and the 
exponential correlation length in the low temperature phase of the Ising spin 
model is equal to the inverse mass of the 0"*"^ glueball in confined phase of 
the Z2 gauge model. While there is no direct experimental particle physics 
relevance of this result, it is interesting for theoretical reasons to compare 
mo++ / obtained from different gauge theories. Finally we also provide an 
updated estimate of the finite temperature transition T^/ y/a. Note that here 
we refer to the temperature of the two-dimensional quantum field theory. Its 
temperature is given by T = l/(aLo) and should not be confused with the 
temperature of the three-dimensional classical system defined above. In the 
following we shall denote the critical value of Lq by A'^; i.e. = \/[aNt). 

The content of this paper is the following: In section [2] we define the 
interface free energy for finite interface area Li x L2 and finite transverse size 
of the system Lq. Next, in section [31 we briefly summarize the predictions 
for the dependence of the interface free energy on (Li,L2), according to 
an effective string-like description. In section [4] we present our numerical 
method to compute the interface free energy. Our results for square and for 
rectangular interfaces are presented in section [5l while section [6] contains our 
results for the universal amplitude ratios. A summary and our conclusions 
are given in section [71 The numerical integration methods are presented in 
the appendix [21 

2 Definition of the interface free energy 

The basic quantity that we shall determine numerically is the ratio between 
the partition functions of the system with anti-periodic and periodic 
boundary conditions Z^. The purpose of this section is to provide a definition 
of the interface free energy in terms of this ratio. 

The ratio Z^/Z^ can be expressed in terms of the eigenvalues A„x of the 




(4) 



3 



transfer matrix and the parity p^^ — ±1 of its eigenstates: 



(5) 



For Lq 3> ^, where ^ = — l/ln(Ais/Aoa) is the bulk correlation length or the 
inverse of the mass of the theory, the partition function ratio in eq. ([5]) is 
dominated by the largest eigenvalues Aqs and Aoa: 

^p Ao^° + A^: 1 + (Aoa/Ao.)^o ■ ^ ^ 

In this regime, the so-called tunneling mass: 

mt = -ln(Aoa/Aos) (7) 

can thus be obtained from: 

1 fi-zjzA 

mt = --f-ln[ ^ „ ,^ . (8) 



Lo \1 + ZJZ, 

Now let us relate the ratio of partition functions with the phenomenological 
picture of interfaces separating the phases of positive and negative magneti- 
sation. We assume that, to the leading approximation, the free energy of an 
interface is proportional to its area. Hence, for finite Lq, in the Li, L2 — > oo 
limit, there is only one interface in the system with anti-periodic bound- 
ary conditions and none in the system with periodic boundary conditions. 
Based on this scenario, the interface free energy is naturally defined as: 

F« = -ln(Z,/Z,)+lnLo , (9) 

where the IuLq term takes into account the "entropy" due to the fact that 
the interface can be located at any point in the xo-directionH 

Note that for finite ^1,-^2 the value of F^^^ depends on Lq and in par- 
ticular, the limit Lq ^ cxd is not finite. This last problem is related to the 



^Eq. ([SJ can be justified as it follows: In tfie basis of slice configurations E, the matrix 
associated with anti-periodic boundary conditions is given by Pe',e = ^e',-E) where — E 
means that all spins in the slice are flipped. Since the external field h is vanishing, the 
transfer matrix commutes with Pe'.e; furthermore, P squares to the identity, therefore 
it has eigenvalues pnx = ±1. We label eigenvectors with p„x = 1 by x=s and those with 
p„x = — 1 by x=a. Following standard conventions Anx is decreasing with increasing n. 

^In principle, one might also add a further In 2 term, to take into account that the 
positive magnetization domain can be realized on the left-hand side of the interface and 
the negative one on its right-hand side, or vice versa. 
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fact tliat for sufficiently large Lq, it is favoured by the entropy to create 
additional pairs of interfaces. 

The presence of additional pairs of interfaces can be addressed in the 
dilute gas approximation. I.e. we assume that the interaction of two inter- 
faces is short ranged and that the average distance between interfaces is large 
compared with the range of the interaction. For n separate, non-interacting 
and indistinguishable interfaces with the free energy one obtains: 

Zi = Y^ exp(-nF,) = ^ exp[-n(F, - IuLq)] . (10) 

n ' n ' 

The sum runs over non-negative even integers in the case of periodic bound- 
ary conditions, and positive odd integers in the case of anti-periodic bound- 
ary conditions, and the ^ factor takes into account that the interfaces are 
indistinguishable. Hence: 

_ Em=o T2^^. exp[-(2m + l)(F, -InLp)] 

Em=o exp[-2m(F, -InLo)] 

= tanh{exp[-(F, - InLo)]} . (11) 

The solution of this equation with respect to F^ provides us with a second 
definition of the interface free energy: 

Upon comparison between eq. (jl2p and eq. dH), the tunneling mass mt can 
be expressed in terms of the interface free energy as: 

= 2exp(-F(2)) , (13) 

which confirms that the definition F^'^\ in contrast to F^^\ has a finite, 
meaningful Lq ^ oo limit. This limit is well approximated for Lq ^ ^. 
All our simulations are done in this regime. Note that for Lg ^ Ct) i-e. 
Za/Zp close to zero, i^J^-* is a good approximation of F^'^\ For most of our 
simulations, this condition is satisfied. 

3 Interfaces in gauge theory: The effective string 
perspective 

In quantum gauge theory, the low-energy behaviour of a confined pair of 
static sources at a distance r might be described by an effective string. In 
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the confining regime, the flux lines between the two sources are squeezed 
into a thin tube, which might be ideahzed as a uni-dimensional object. 
The long-distance properties of the system are dominated by the transverse 
fluctuations of this tube; in this regime, the excitation spectra of the fields 
in the interior of the tube are expected to be much higher-lying. 

Under this assumption, the properties of the system are described through 
a string partition function, obtained integrating over the possible world-sheet 
configurations. Each of them has the topology of a cylinder, and contributes 
a Boltzmann-like factor, whose exponent is given by an effective string ac- 
tion. 

In principle, the functional form of the latter is unknown, however it can 
be constrained, by requiring that it satisfies certain self-consistency proper- 
ties, and that it yields the correct physical limit for large distances r between 
the two sources. 

This approach underlies the models that have been proposed by Polchin- 
ski and Strominger in the 1990's [16]: 



\id+x-d-X) 



D - 26\ (dlx ■ d-X)(d+X ■ dlx) ^, 
' h 0\r 



12 ) • d-Xf 



(14) 



(in which are light-cone world-sheet coordinates, and a is a length scale 
related to the string tension) and by Liischer, Symanzik and Weisz already 
at the beginning of the 1980's [17, 18] (see also [19, 20] for more recent 
developments) : 

S'efi = arL /iL 5o + Si + 52 + . . . , with: So = ^ j d^i{dahdah) , 

(15) 

(where L denotes the length of the closed world- lines of the static sources, 
a is the string tension, and /i is a coefficient associated to a perimeter-like 
term). The construction of a string action in the form of eq. (jl4p allows a 
generic conformally invariant world-sheet QFT, with the coefficient of the 
various terms fixed by anomaly cancellation in any dimension D. The action 
can be built by converting the path integral for the collective coordinates of 
the underlying field theory to covariant form. The X field is unconstrained 
(and the model still represents a generic interface in D dimensions), but the 
term appearing in the second line of eq. ()14p takes the Polyakov determinant 
in conformal gauge into account. Eq. (jl4p yields an effective string model 
which can be expanded around the long-string vacuum. Poincare invariance 
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constrains the 0{r^'^) term of the string spectrum to the form it has in the 
Nambu-Goto model [21-24] (for the definition of the Nambu-Goto model 
see below). 

In eq. p^ . the 5*0 term describes a conformal model, while the other, 
higher-dimensional, Sn terms are responsible for the string self-interactions. 
In this case, the variable h represents a vector with D — 2 components, 
that describes the fluctuations transverse to the reference plane (physical 
gauge). It is interesting to note that the first term beyond So, which is 
expected to be a "boundary term" [19] , is actually forbidden by open-closed 
string duality [20] . 

Among the main implications of the effective string model, we mention 
the existence of a negative, 0(r~^) correction to the asymptotic linear po- 
tential (the Liischer term, which is a Casimir effect), and the logarithmic 
growth of the square width of the flux tube [25-27]. Both aspects are re- 
lated to the fact that, at leading order, the effective string fluctuations in a 
L'-dimensional target space can be modelled as {D — 2) free, massless bosons. 

Since the infinite number of terms appearing on the right-hand sides of 
eq. (I14|) and eq. (llSp are not known a priori, it is, in general, not possible 
to work out all-order predictions for the observables of interest. 

Alternatively, one might use an explicit ansatz on the functional form of 
the effective string action (consistent with the constraints mentioned above, 
and compatible with the other effective models at its lowest orders): this 
allows to address the complete mathematical calculation for the expectation 
values of the physical observables, and to perform an all-order comparison 
with the numerical results. 

A natural choice for the effective string action (for any world-sheet ge- 
ometry) is the area of the string world-sheet itself. For the case of a closed 
interface, it can be expressed introducing the ^ coordinates over the inter- 
face, and the ga/3 metric induced by the embedding in the target space: 



(the string tension a has energy dimension 2). 

Eq. (jl6p has a natural interpretation in the context of string theory, 
where it represents the (Euclidean space formulation of) the model due 
to Nambu and Goto [28,29], describing the relativistic quantum dynamics 
of a purely bosonic string. Although it is well-known that this model is 
affected by an anomaly (breakdown of rotational symmetry out of the critical 
space-time dimension D = 26) and is non-renormalizable (because it is non- 
polynomial), it has been studied as a possible effective description of the 




(16) 
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low-energy dynamics of confining gauge theories. The reason is that, in 
the infra-red regime, the lowest-lying degrees of freedom associated with 
a confining flux tube are transverse fluctuations, and are expected to be 
modeled by a bosonic string-like dynamics; in that case, a represents the 
asymptotic value of the string tension. In particular, the geometry of an 
interface with periodic boundary conditions in both directions would be 
associated to the description of a torelon, i.e. a string winding around a 
compact target space, that has been studied in ref. [30]. 

A number of implications for this effective description of confinement 
have been derived theoretically and tested numerically in the literature [31- 
51]. 

On the other hand, in a condensed matter physics context, eq. p6]) 
corresponds to the "capillary wave model" [1]; at first order, it describes the 
transverse fluctuations of a membrane as free, independent, massless modes, 
whereas the subleading terms introduce (self-) interactions. In this context, 
a can be interpreted as an (asymptotic) interface tension, which does not 
depend on the local orientation of the normal to the inflnitesimal surface 
element. 

A perturbative expansion in powers of {aLiL2)~^ yields the following 
result for the partition function of the interface with periodic boundary 
conditions in both directions [44-46]: 



A -2 
Z = — exp(-crLiL2) V (iu) /r] [i] 
'u 



(TL1L2 



1 



{'JL1L2Y 



(17) 

where the parameter A can be predicted invoking a perturbative argument 
for the (j)'^ scalar fleld theory in three dimensions, t = iu = iL2/Li is the 
modulus of the torus associated with the cross-section of the system, ry is 
Dedekind's function: 



7/(r) 



1/24 



and f{u) is defined as: 



11(1-9") > (? = exp(2^^r) , 



n=l 



UE2 {iu) 



vr 3 
-UE2 (^u) + - 



where E2{t) is the first Eisenstein series: 



(18) 



(19) 



E. 



:(r) = l-24j]- 



n q 



n=l 



q = exp(27rir) , 



(20) 



8 



In particular, for u = L2/L1 = 1 one gets /(I) = 1/4. 

The interface free energy for square lattices of size Li = L2 = L takes 
the form: 



F=aL'^-lnX-—^ + 0{—^] . (21 



1 „ / 1 



4aL2 V(f^^^)^ 
This is the theoretical expectation which, in the following section, will be 
compared with our numerical results for — see eq. ()12p . 

More recently, a different approach to calculate the partition function 
was proposed in ref. [15]: this method is more elegant and powerful with 
respect to the perturbative expansion in powers of {(jLiL2)^^ , and it takes 
advantage of the standard covariant quantization techniques for the bosonic 
string. The power of this method relies in the fact that it allows to resum 
the complete loop expansion for the interface partition function at all orders, 
with a final result for the interface partition function in d dimensions I^"^^ 
taking the form of a series of Bessel functions: 



I^''^ = 2(^—y FrV^ J^^CfeC^-fc (-) K^{aA£) , (22) 



m=0 k=0 



(where A = L1L2, and Vr is the product of the system sizes in the transverse 
directions) and a consistent, closed-form expression for the spectrum levels: 



£ = £,,^ = ^l + ^^(m-i^y^^(2k-,nf , (23) 

that agree with those presented in ref. [41]. 

For the case of an interface with the boundary conditions of a cylinder, 
an analogous result was derived in ref. [52], while the associated energy 
spectrum had already been known since the Eighties [53]. 

In eq. (j23p . for d = 3 and m = k = 0, the argument of the square 
root becomes negative for aLf < 7r/3. This is known as the tachyonic 
singularity in the effective string framework (see [54]) and can be physically 
interpreted as the signature of a high temperature deconfinement transition: 
For temperatures higher than T^/ \fa = \/2>/t^ the string vanishes and quarks 
are no longer confined. Equivalently, in the dual model, the Ising spin 
model, for Li < Nt the interface tension vanishes and a transition from a 
ferromagnetic to a paramagnetic phase occurs. 

This interpretation was discussed for the first time by Olesen in the 
case of Polyakov loop correlators [54] and holds essentially unchanged in 
the present case, although the Arvis spectrum [53] on which Olesen's result 
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was based is very different from the one we have in the present case — see 
eq. (|23p . In fact, the lowest state (the one which drives the phase transition) 
is the same in both spectra. 

Notice that one should not expect that this analysis provides an exact 
result for the finite temperature transition. Indeed the Svetitsky-Yaffe con- 
jecture [55] suggests that the finite temperature transition in the Z2 gauge 
theory belongs to the universality class of the 2D Ising spin model. On 
the other hand, the analysis of the string picture gives mean-field critical 
exponents (i.e. v = 1/2). Once more this observation indicates that the 
Nambu-Goto effective string should be better considered as a mean field 
description, which is particularly effective at low temperatures and/or large 
distances. 



4 Method to compute the interface free energy 

Let us define the interface energy as: 

E^ = E^-E, , (24) 

where (Ep) is the expectation value for the energy of a system with 
anti-periodic (respectively, periodic) boundary conditions. 

The internal energy of a system is given by the derivative of the reduced 
free energy with respect to /3: 

dlnZjP) _ E{s}^M-m{{s})] H{{s}) 

dp E{s}^M-PH{{s})] ■ ^ > 



From eq. (j25]l . by integration from (3q to (3 it follows: 
By adding InLo on both sides of the equation we get: 



In — 



= - In — 

/3 -^P 



+ / d/3 EXP) . (26) 
J0n 



FP{P)=FPiPo)+ ["dp E^P) . (27) 



A) 

In general it is difficult to determine free energies directly in a single 
Monte Carlo simulation. On the other hand, expectation values such as Ep 
and E^ can be easily determined. 

It is rather an old idea to compute free energies (in particular: interface 
free energies) from eq. (|27p — see, e.g., the first few references in [56]. 



10 



For /?o there are different possible choices: For (3 < (3^ the interface 
tension vanishes, and with a suitable choice of /3o < /3c (depending on the 
interface area) F^^\I3q) vanishes to a very good approximation. Alterna- 
tively, one might start the integration from large values of /3 ^ /3c, where 
F^^\Pq) can be obtained from the low temperature expansion. 

Here we follow the strategy discussed in [56]: F/^)(/3o) is computed using 
the boundary flip algorithm [57] at a /3o value corresponding to a Z^/Z^ 
ratio of the order of 1/10. For such a choice, FJ^)(/3o) can be accurately 
determined using a moderate amount of CPU time. 

In practice, we have performed simulations for a finite number of inverse 
temperatures (3o < f5i < P to obtain values for E^{j3i). The integration (p7|) 
must then be performed by using some numerical integration scheme. For a 
detailed discussion of the schemes that we have used see the appendix El 

In principle, the numerical integration, along with its (small) systematic 
error, could be avoided. One might compute: 

F(/? + A/?) - F(/?) = -log[Z(/? + A/?)/Z(/3)] 

= -log{(exp[-A/?i/({4)])^} (28) 

in a Monte Carlo simulation that generates the Boltzmann distribution cor- 
responding to (3. Here A/3 has to be chosen such that A/3-^ {H"^) " (H)"^ is 
of order one, to keep the statistical error under control. Then -F(/3) — F{Pq) 
is computed by a sequence of such A/3 steps. 

For an interesting new alternative, using a generalized Jarzynski relation, 
see ref. [58]. 

Quite a different strategy to compute is based on the so-called snake 
algorithm [59,60]: A sequence of boundary conditions is introduced, that 
interpolates between periodic and anti-periodic boundary conditions. The 
boundary in the 0-direction is progressively filled with J(xy) = ~1- I^i 
refs. [3-14] we have used this approach to compute ratios of Polyakov-loop 
correlators. 

In the present work, however, we do not use the snake algorithm, since 
for boundary conditions close to the anti-periodic boundary conditions the 
simulations are very difficult. If the last layer of bonds is not completely 
filled with antiferromagnetic couplings J{xy) = — 1 yet, then it is energetically 
favorable that the interface sticks to the boundary. On the other hand, there 
is an entropy gain, when the interface moves freely along the 0-direction. 
At some stage these two effects are approximately in balance, resulting in 
extremely long auto-correlation times for any known algorithm. 

A major reason for using the numerical integration is that, given the 
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large number of lattice sizes and values of /?, it provides the most efficient 
way to organize the simulations and to keep the resulting data under control. 

4.1 Monte Carlo simulations 

In order to compute the starting-point free energy F^IPq), we have used a 
variant of the boundary flip algorithm [57] — see refs. [13, 61] for a dis- 
cussion. The update is performed using the single cluster algorithm [62]. 
Typically, we have performed O(IO^) up to O(IO^) measurements. For each 
measurement, mostly 10 single cluster updates were performed. 

In order to compute the energy for the systems with periodic and anti- 
periodic boundary conditions, we have used a demonized local Metropolis 
algorithm implemented in multispin coding technique: details can be foTind 
in ref. [63]. This way, n^it systems run in parallel in our implementation 
(^bit = 32 or ribit = 64, depending on the machine used). Most simulations 
were performed with the local algorithm alone. 

For each measurement, 12 complete update sweeps were performed. For 
n^it = 32 we performed either 100,000 or 96,000 measurements for each copy 
of the system and each value of /3; for n^.^ = 64, 50,000 measurements for 
each copy of the system were performed. In order to ensure equilibration, we 
have taken 24,000 local update sweeps. Note that this is an overkill for the 
smaller lattice sizes and the larger /? values. However, due to the enormous 
number of individual simulations, we could not check carefully each of the 
runs. Therefore we decided to use a common number of thermalization 
updates, that is suitable for all of the parameter choices we considered — 
including the most difficult cases. 

For P values close to Pq, there is a non- negligible probability that more 
than (periodic boundary conditions) or 1 (anti-periodic boundary condi- 
tions) interfaces are formed in the system. Most likely, the local update 
is not capable of generating these additional interfaces within the given 
number of update cycles. Therefore, we performed single cluster updates 
in addition. In this step, we could not make use of the multispin coding 
technique, therefore the cluster update was performed one-by-one for the 
ribit systems. One single cluster update is performed per measurement; we 
made no attempt to optimize the ratio of cluster and local updates. For a 
discussion about a similar combination of algorithms, see ref. [64]. We have 
used this method, instead of the local update only, for (3 = f3o + mAp with 
m < 10. 

Going to larger interface areas than those studied in the present work, 
it would be advisable to use the interface cluster algorithm introduced in 



12 



ref. [65] and further discussed in ref. [66]. 

In table [1] we give a summary of our runs for square interfaces and in 
table [2] for the asymmetric interfaces. In total, order of thousands individual 
simulations were performed. 

4.2 Numerical integration 

The numerical evaluation of integral (j27p was done using standard numerical 
integration schemes which are summarized in the appendix El 

All the schemes that we have considered can be written in the form: 

N 

FP{P) = FW(p,) + Y^c,Af3EM+jm+0{N-^') , (29) 

j=0 

where A/3 = {(i — (3o)/N and Xlj^o'^i ~ ^^'^ final estimates, we 

have used schemes with an 0{N~'^) integration error. In order to get a 
quantitative estimate of the integration error, we have compared the re- 
sult obtained from different schemes; e.g. scheme ()A.3P and scheme ()A.4p . 
Furthermore, we have performed the numerical integration for the theo- 
retical predictions of F^{Li, L2, a), as discussed in section [3l along with 
a{P) = (To(/3 — Pc) X [1 + b{P — + c{(3 — /3J], with coefficients similar 
to those reported below. We found that the error of the integration is at 
least two orders of magnitude smaller than our statistical error, and is hence 
ignored in the further analysis of the data. 

4.3 Propagation of the statistical error 

The statistical error e of F^^^ (/?) is computed using standard error propaga- 
tion: 

e^F^m = e2[F(i)(/3o)] + (A/?)^ c|e2[i?,(/3o + jA/3)] , (30) 

i 

where e is the statistical error and the Cj coefficient is given by the integration 
rule. 

In order to get correct fits for F^ at different values of /?, we have to 
evaluate the covariances of F^ at different values of (3. Let us consider /?2 > 
/3i: Due to the fact that fP{I3q) and E,{P) with P < Pi are obtained in a 
common set of simulations, F^^\Pi) and F^^\P2) are statistically correlated. 

The covariance is defined as: 

coy{A,B) := {[A - {A)][B - {B)]) . (31) 
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6.55408(116) 
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6.77107(123) 
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Table 1: Summary of the simulations for the square interfaces. For each 
linear extension L = Li = L2 of the interface, the table gives the starting 
point Po of the integration, the maximal inverse temperature /J^ax that has 
been simulated, and the step-size A/3. In the case of L = 17, 18, 19, 20 we 
have two intervals with different A/3. The initial value of the integration 
F^^^ (Po) has been computed with the boundary flip algorithm. For details 









24 


0.222402 


5.19305(69) 


28 


0.222402 


5.53093(93) 


32 


0.222402 


5.87624(96) 


36 


0.222402 


6.21529(116) 


40 


0.222402 


6.54949(141) 


44 


0.222402 


6.87239(137) 


48 


0.222402 


7.19102(163) 



Table 2: Summary of runs with asymmetric lattices Li ^ L2, in the same 
notation as in the previous table. Lq = 96, L2 = 64, A(3 = 0.00005 and 
/3^ax = 0.226102 throughout. 

In our case: 

fP{P) = F('\Po) + ^c,{P)E, , (32) 
j 

where the Ej and F^^\Pq) are statistically independent. Hence: 

cov(F(i)(/3i),Fii)(/32)) = var(F(i)(/3o))+^c,(/3i)c,(/32)var(£;,) 

j 

« vaiiPPiPi)) . (33) 

The last equality is only approximate, due to the fact that Cj(/?i) 7^ 1 for 
the last few j < m. In the limit A/? — > 0, the approximation becomes exact. 
In our data, we have checked that the approximation is very good, and it is 
therefore used in the fits. 

5 Numerical results for the interfaces 

5.1 Square interfaces 

First we have analysed the data for the square lattices. In comparison 
to our previous work ref. [13], we have results for more than four times 
larger interface areas. This allows us to use the interface tension as a fit 
parameter, while in ref. [13] we had to take it from ref. [7], where Polyakov- 
loop correlators were studied. 

Furthermore, here we have data for a large range of inverse temperatures 
allowing us to address the question of corrections to scaling. 
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As a starting point, let us first discuss the results for /? = 0.226102. Note 
that for (3 = 0.226102, the finite temperature phase transition occurs at Lq = 
8 [67]. In table[3]we have summarized our results for the interface free energy 
at /? = 0.226102. Note that here we have converted, using eqs. (j9ll2p . 
F^^\ which is the result of our numerical integration, to F^'^\ which is less 
dependent on Lq and is therefore more suitable for the comparison with the 
theoretical predictions. For Lq/L fixed, the difference between F^"^^ and F^^^ 
goes down exponentially as the interface area increases. For our numerical 
results at P = 0.226102 this difference is larger than the statistical error 
only for L < 20. 

Similarly to ref. [13] we have fitted the data with the ansatze: 

fP =aL' + co + ^ (34) 

and 

^.)=,.^ + e„ + ^ + ^. (35) 

where a, cq, C2 and C4 are the free parameters of the fits. At this stage of 
the analysis we made no attempt to compare with the full NG-prediction 
which can be obtained from eq. (j22p . 

Results of fits with the ansatz (|34p are given in table [H Starting from 
-^min = 19, the x^/d.o.f. is smaller than one. The fit result for C2 is, up 
to Lniin = 24, decreasing with increasing L^i„. For L^i^ = 24 we get C2 = 
—0.246(13) which is consistent with the NG prediction C2 = —0.25. Since 
the fit result is increasing with L„,in we might consider C2 = —0.246(13) as 
an upper bound. 

As a check we have performed fits with the ansatz (I35p : the results are 
summarized in table [5l The x^/d.o.f. is below one for all L^in available. 
C2 is now increasing with increasing L^^^. Unfortunately, no stable estimate 
for C4 is obtained. Higher order corrections seem to play an important 
role. The results C2 = —0.280(38) from L^^^ = 20 might serve as lower 
bound for 02- Combining the results of the fits with the ansatze (I34|35p we 
might summarize our results as: —0.246(13) > C2 > —0.280(38), which is 
fully consistent with the theoretical prediction. As our final result for the 
interface tension we take (t(0. 226102) = 0.0105255(11), obtained from the 
ansatz (I34p and L^in = 24. The comparison with results from the ansatz (|34p 
suggests that systematic errors should not be larger than the statistical error 
that is quoted. 

We have repeated this type of analysis for f3 = 0.223102, 0.223452, 
0.223952, 0.224752, 0.227202, 0.228802, 0.230002, 0.236025, 0.24 and 0.24607. 
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L 


Fj2) 


18 


6.00956(34) 


19 


6.40442(38) 


20 


6.81999(45) 


21 


7.25617(50) 


22 


7.71334(51) 


23 


8.19024(56) 


24 


8.68757(58) 


25 


9.20809(77) 


26 


9.74659(78) 


27 


10.30706(84) 


28 


10.88919(87) 


29 


11.48975(92) 


30 


12.11320(98) 


31 


12.7558(10) 


32 


13.4210(11) 


33 


14.1074(12) 


34 


14.8145(13) 


35 


15.5415(13) 


36 


16.2881(15) 


38 


17.8477(16) 


40 


19.4919(17) 


44 


23.0292(22) 


48 


26.9095(24) 


52 


31.1193(28) 


56 


35.6618(32) 


64 


45.7769(40) 



Table 3: Interface tension F^^ for square interfaces at /3 = 0.226102 as a 
function of the linear lattice size L. 
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18 


0.0105283(7) 


2.6589(11) 


-0.209(4) 


1.40 


19 


0.0105273(8) 


2.6611(13) 


-0.219(5) 


0.95 


20 


0.0105267(8) 


2.6625(15) 


-0.226(6) 


0.81 


21 


0.0105264(9) 


2.6635(17) 


-0.230(7) 


0.79 


22 


0.0105262(9) 


2.6640(19) 


-0.234(9) 


0.82 


23 


0.0105257(10) 


2.6654(22) 


-0.242(11) 


0.76 


24 


0.0105255(11) 


2.6661(25) 


-0.246(13) 


0.78 


25 


0.0105262(12) 


2.6636(30) 


-0.229(17) 


0.69 


26 


0.0105259(13) 


2.6649(34) 


-0.239(21) 


0.70 



Table 4: Fits with ansatz (I34p for square interfaces at /? = 0.226102. 

Throughout we find for fits with the ansatz (I34p that the numerical result 
for C2 is decreasing with increasing -Lmin; until it starts to fluctuate. In the 
case of the fits with ansatz (]35p we see that C2 is increasing with increasing 
L^i„ for P < 0.230002; for larger values of P it decreases. For /? < 0.230002, 
our results are consistent with the theoretical prediction (for the scaling 
limit) C2 = —0.25. For larger values of /?, deviations become visible: for 
instance, for /? = 0.24, the result from a fit with ansatz ([3^ and L^i„ = 12 
is C2 = -0.31(2). 

In order to disentangle corrections to scaling from truncation effects in 
the ansatz, we have studied C2 as obtained from fits with the ansatz (I34p 
and \/o"-L„i„ 2 fixed. Naively fitting all data for /3 < 0.24 we get: 

C2|l^,„=2/V? = -0.217(4) -0.98(15)a ' (^6) 

with x^/d.o.f.= 0.55. In figured] we show our data along with the result of 
this fit. We conclude that C2 is affected by corrections to scaling. However 
the corrections are, within the numerical precision, proportional to o" oc 
i.e. they vanish much faster than where uj = 0.821(5) [68] is the 

exponent of the leading correction to scalingH This could be explained by 
the fact that the effective interface model only assumes that the symmetries 
of the continuous space are restored. The restoration of rotational symmetry 
is indeed associated with a correction exponent oj' « 2 [71]. 

■^Field theoretical methods give slightly smaller values: lo = 0.814(18) from the e- 
expansion and lu = 0.799(11) from perturbation theory in three dimensions fixed [69]. In 
ref. [70] the value lj = 0.845(10) was obtained from Monte Carlo simulations of the <f>* 
model on the lattice. 
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Figure 1: Results for C2 from fits with the ansatz (I34p and L^in = 2/\/a. 
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18 


0.0105246(11) 


2.6701(29) 


-0.299(22) 


0.20(5) 


0.71 


19 


0.0105250(12) 


2.6689(35) 


-0.288(29) 


0.17(7) 


0.72 


20 


0.0105252(14) 


2.6682(42) 


-0.280(38) 


0.14(10) 


0.75 


21 


0.0105252(15) 


2.6678(51) 


-0.275(50) 


0.13(15) 


0.79 


22 


0.0105252(15) 


2.6679(56) 


-0.276(60) 


0.14(19) 


0.84 



Table 5: Fits with ansatz (I35p for square interfaces at /? = 0.226102. 
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Figure 2: The constant Cq as a function of (3. The quantity is defined in 
eq. (j37|) in the text. The value for P = 0.27604 is taken from ref. [13]. 

Next we study the behaviour of cq. In the scahng hmit, this quantity 
should behave like: 

co{P) = Co-\ln[am. (37) 

In fig. [2] we have plotted our results for Cq as a function of (3. 

Within our numerical precision there is no sign of corrections to scaling 
whatsoever for /? < 0.24607. Only for P = 0.27604, the value is taken from 
table 5 of ref. [13], a clear deviation is visible. Unfortunately, we have no 
clear theoretical understanding why corrections to scaling should be so small 
in this quantity. 

Results for the interface tension are summarized in table [H All these 
results are taken from fits with the ansatz (j34p using L^i,^ ~ 2.5/y^. As 
discussed above, for the case (3 = 0.221602, systematic errors should be 
smaller than the statistical error that is quoted. 
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5.2 Global fit of the data for the interface free energy 

In the neighbourhood of the transition, the interface tension behaves as: 

a{(3) = aof" X (1 + + bt + ct^^ + df + ...) , (38) 

where t = (3 — (5^ \s the reduced temperature. The most accurate result 
for the inverse critical temperature is f3^ = 0.22165455(5) [68]. The critical 
exponent /U of the interface tension is related with the critical exponent of 
the correlation length as /U = 2u. The most accurate values given in the 
literature are v = 0.63012(16) [72] from the analysis of high temperature 
series expansions, and v = 0.63020(12) [68] from a finite size scaling analysis 
of Monte Carlo data. The same study provides oj = 0.821(5), hence 6 = 
uio = 0.5174(33). Note that the value of 6' = 1.05(7) [73] has quite a large 
uncertainty. 

Since 20 ~ 0' ~ 1, we take as ansatz for the fits: 

a{p) = aof" X (1 + at'^ + bt) . (39) 

Fitting our new results for cr(/3) to this ansatz would require to take into 
account the cross-correlations among the values of (t{P) for different values 
of p. Instead of computing these cross-correlations, we performed fits for 
fitting the L and the /3 dependence at the same time. The cross- 
correlations of at different values of /? can be easily obtained as discussed 
in subsection 14.31 

Based on the results obtained above, we performed a four parameter fit 
of the data for F^{L,f3). To this end, we have used the ansatz: 

F.(L,/3) = + Co -^Ham - , (40) 

where the interface tension is given by the ansatz ()39p . namely, the free 
parameters are do, a, b and Cq. The critical exponents and the inverse 
critical temperature are fixed by their best estimates given in the literature, 
which are quoted above eq. ([39|) . 

In the fit, we have used results for the interface free energy at the 
same values of j3 as discussed in the previous subsection: j3 = 0.223102, 
0.223452, 0.223952, 0.224752, 0.227202, 0.228802, 0.230002, 0.236025, 0.24 
and 0.24607. Our data would allow to use more values of (3. However, little 
information would be added this way, since, by construction, the interface 
free energies at close-by values of [3 are highly correlated. 

After some experimenting we decided to take our final estimate from 
a fit with input data characterized by /? < 0.227202, and F^ — \nL > 8, 
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which roughly corresponds to \faL > 3. In total, 51 data-points satisfy 
this criterion. The results for the fit parameters are ctq = 10.083(8), a = 
-0.479(26), b = -2.12(19) and Co = 0.3895(8), where xVd-o.f. = 0.79. 

In order to check the L dependence of our result, we have repeated the 
fit with P < 0.227202, and — InL > 4 (which corresponds roughly to 
y/aL > 2.1) and L < L,„ax = 44. This means that the range in L is roughly 
\/2 times smaller than that of the previous fit. In total, 84 data points satisfy 
this criterion. The results of this fit are do = 10.080(8), a = -0.471(24), 
b = -2.20(17) and Co = 0.3915(4) with xVd-o.f. = 0.89. 

Next, we changed the /3- interval of our fit. We included data that satisfy 
the criteria 0.224302 < /? < 0.233, - InL > 8, L < L^,, = 44. There 
are 55 data points that satisfy these criteria. The results of the fit are 
ao = 10.085(5), a = -0.484(12), b = -2.08(7) and Co = 0.3886(6) with 
xVd.o.f. = 0.64. 

The results for ao are consistent among the three different fits. The 
differences, possible within the statistical error, of these results provide an 
estimate of the systematic error due to finite-L and large-t corrections that 
are not taken into account by the ansatz. We arrive at: 

(70 = 10.083(8) [26] + (41) 
22330(/3, - 0.22165455) + 174(i/ - 0.6302) - 0.237(6* - 0.5174) 

and 

a = -0.479(26) [120] (42) 
-55866(/3, - 0.22165455) - 149(z/ - 0.6302) - 2.82(0 - 0.5174) . 

The number in the brackets gives the systematic error caused by corrections 
to the ansatz, as discussed above. In the second line of eqs. (I41I42P we give 
the dependence of the result on the input parameters /3c, and 9. The 
dependence of Co on (3^, v and 9 is small enough to be neglected. We take: 

Co = 0.3895(8) (43) 

as our final result. The comparison of the three fits done above suggests that 
the systematic should not be larger than the statistical error. In table [6] we 
compare results for a obtained from the global fit with the results obtained 
from analysing single values of /? in the previous subsection. For all values 
of 13 the results are consistent. 
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/9 
P 


(T global fit 


(T 


n OOQ1 no 


U.UUzDUoo(bj( / ) 


U.UUzdU4o(0o j 




U.UUo4i io[o)[o) 


U.UUo4i0z(oi j 




n nn4f!'^Q7i'fiVi i 

u.uu^uoy / i^u j \^ 1 1 J 




0.224752 


0.0067258(6)(16) 


0.0067269(17) 


0.226102 


0.0105254(7)(28) 


0.0105255(11) 


0.227202 


0.0138217(8)(42) 


0.0138220(17) 


0.228802 




0.0188659(13) 


0.230002 




0.0228068(12) 


0.233 




0.033114(15) 


0.236025 




0.044019(9) 


0.24 




0.058913(5) 


0.24607 




0.082510(5) 



Table 6: Final results for a at given values of the inverse temperature 
/5. In the second column we give a as obtained from our global fit. The 
statistical error of a is properly computed, and the second error quoted is 
the systematic one. It is estimated from comparing results from different fit 
ranges. The third column gives a obtained from fits with the ansatz (j34p 
and Lmin 2.5/v^. For the three smallest values of (3 the global fit provides 
more accurate results for a than the fits with the ansatz ()34p . 

5.3 Comparison with the hterature 

Using the definitions i = {(3 — (3^) / 13^ and a = a^t^ x (1 + ...) we get: 

ao = 1.510(4) , (44) 

where the error is dominated by the uncertainty of v. Note that this result 
is perfectly consistent with (and more precise than) the most accurate result 
(To = 1.50(1) given in the literature [74] using Monte Carlo data of ref. [63]. 
A more comprehensive list of results for is given in table 8 of ref. [66] . 

Also our results for the interface tension at given values of [3 can be 
compared with values given in the literature. Here we only give a small 
selection of the most recent results. For more see e.g. ref. [75]. In [7], study- 
ing Polyakov-loop correlators, we find a = 0.0105241(15) and 0.044023(3) 
for [3 = 0.226102 and 0.236025, respectively. These values are completely 
consistent with the estimates of the present work. One should note that 
they are obtained using a completely different numerical procedure. 

In Fig. 4 of ref. [36] the results o- = 0.004782(6), 0.01011(10), 0.022798(2) 
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and 0.02752(10) for /3 = 0.224, 0.226, 0.23 and 0.23142 are provided. 
These can be compared with a = 0.004761(2), 0.010228(4), 0.022800(1) and 
0.027603(2) for the same values of /?, taken from our global fit. Our values 
are consistent with those of ref. [36], except for /3 = 0.224, where we observe 
a discrepancy by three and a half standard deviations. 

In ref. [66] the interface tension has been computed in a similar way as 
in the present work. Our present results are by approximately a factor of 
ten more precise than those of ref. [66]. The results quoted in ref. [66] are 
consistent with our present estimates within two standard deviations. 

The results of ref. [44] , using the two- loop approximation to fit the data, 
are a = 0.004778(14), 0.006547(69), 0.009418(61) and 0.014728(40) for /3 = 
0.2240, 0.2246, 0.2258 and 0.2275 to be compared with a = 0.004761(2), 
0.006319(2), 0.009418(61) and 0.014740(5). While the results for f3 = 0.2240 
and 0.2275 are perfectly consistent, there is a mismatch by 3.3 and 3.8 
standard deviations in the case of /3 = 0.2246 and 0.2258, respectively. This 
is likely due to the fact that for these two values of /? only small interface 
areas were available and too small areas had been included into the fit. 

5.4 Square interfaces: small L 

In figure Elweplot Fp^-aL'^+0.5ln{a) function of \faL. The numerical 
values for a are taken from table [6l Note that these values of a are obtained 
from rather large values of \faL (i.e. are little affected by higher order 
corrections). In addition to the numerical data for square interfaces at /3 = 
0.223102 and /? = 0.226102 (roughly corresponding to the critical values 
of /9 for which the finite temperature transition occurs for A^^ = 16 and 
Nt = 8, where Nt denotes the number of lattice spacings in the "inverse 
temperature" compactified direction) we give the 2-loop prediction and the 
full Nambu-Goto result. In the case of the string predictions we have taken 
Co = 0.3895 into account. We observe that, numerically, for ^/aL ^1.6 
there is very little difference between the two-loop approximation and the 
full NG result. By expanding eq. ()22p we find that the coefficient of the 
1/((tL^)^ term for the full NG result is approximately equal to —0.017. 
Within the statistical error, the Monte Carlo results for the two values of 
P fall on top of the 2-loop and full NG predictions for ^/aL ^ 2.2; note 
that the 1-loop approximation predicts F^ — aL?' to be constant. For smaller 
y^L, the data rather abruptly depart from the string prediction. This 
indicates that (mainly) not 0(1/((tL^)^) but rather higher order corrections 
are responsible for the deviation. Still down to ^/aL « 1.8 the Monte Carlo 
data for the two values of (3 fall on top of each other within the error- 
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Figure 3: Comparison of the 2-loop prediction, the full Nambu-Goto result 
with the data for square interfaces at P = 0.223102 and P = 0.226102. 

bars. For smaller y/o'L, differences become visible, indicating corrections 
to scaling. These scales should be compared with the scale of the finite 
temperature transition ^/aNt ~ 0.81 (for a discussion of this number see 
section El below; the effective string model gives \/aNt = yvr/S « 1.023). 

5.5 Asymmetric interfaces 

In this subsection we compare our results for rectangular interfaces with 
the effective string predictions discussed in sect. [3l Note that the effec- 
tive string results have quite a non-trivial dependence on the aspect ratio 
u = L2/L1. Therefore this comparison is rather a stringent test of the theo- 
retical predictions. We have simulated lattices with the linear sizes Lq = 96 
and L2 = 64. In the remaining direction, the linear size assumes the values 
Li = 24, 28, 32, 36, 40 or 48. We use the values of a and Co obtained in the 
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L2 


Ls 






F«^(/?), full 


(/3), 1-1. 


F:^{P), 2-1. 


24 


64 


96 


6.8887(20) 


6.8855(20) 


6.7495 


6.9974 


6.8347 


28 


64 


96 


7.6935(21) 


7.6929(21) 


7.6537 


7.7875 


7.6821 


32 


64 


96 


8.4626(20) 


8.4625(20) 


8.4518 


8.5380 


8.4632 


36 


64 


96 


9.1996(21) 




9.2012 


9.2632 


9.2062 


40 


64 


96 


9.9227(23) 




9.9231 


9.9713 


9.9253 


44 


64 


96 


10.6203(23) 




10.6278 


10.6674 


10.6288 


48 


64 


96 


11.3138(25) 




11.3209 


11.3548 


11.3213 



Table 7: Free energy of rectangular interfaces, at /5 = 0.223102. The value of 
i?(2)num(-^^ is reported only when different with respect to F/^)""'"(/3). The 
sixth, seventh and eighth column display, respectively, the all- loop (full), 
one-loop (1-1.) and two-loop (2-1.) prediction of the Nambu-Goto string 
model. 





L2 


Ls 






F^'^iP), fun 


F:''if3), 1-1. 


(/?), 2-1. 


24 


64 


96 


7.9988(22) 


7.9985(22) 


7.9397 


8.1043 


7.9802 


28 


64 


96 


9.0175(23) 




9.0066 


9.1016 


9.0212 


32 


64 


96 


9.9965(22) 




9.9962 


10.0594 


10.0023 


36 


64 


96 


10.9407(23) 




10.9455 


10.9917 


10.9482 


40 


64 


96 


11.8682(25) 




11.8707 


11.9070 


11.8719 


44 


64 


96 


12.7728(25) 




12.7803 


12.8102 


12.7808 


48 


64 


96 


13.6734(27) 




13.6791 


13.7049 


13.6793 



Table 8: Same as in table [71 but at /? = 0.223452. 



section 15.11 as input for the theoretical predictions. Therefore we have no 
free parameters to fit, and we can perform a direct comparison among sim- 
ulations and theoretical predictions. We have computed the interface free 
energy Fp^ at /? = 0.223102, 0.223452, 0.223952 and 0.226102, correspond- 
ing to Nt = 16, 14, 12 and 8. Note that only for /3 = 0.223102 and Li = 24 
the difference between F^^ and Fj^^ is slightly larger than the statistical 
error. Our numerical results, along with the various string predictions are 
summarized in tables [7] to [TOl and shown in the plots H] and [5j 

Since we have fixed values of Lq, Li and L2, our data do not allow for a 
check of scaling corrections. Instead, we assume that they are of similar size 
as in the case of the square interfaces and therefore are small for the values 
of P that we consider here. 
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L2 


Ls 




F*(/3), full 


FtiP), 1-1. 


F*(/3), 2-1. 


24 


64 


96 


9.7374(23) 


9.7111 


9.8216 


9.7302 


28 


64 


96 


11.0702(24) 


11.0653 


11.1318 


11.0725 


32 


64 


96 


12.3602(24) 


12.3572 


12.4024 


12.3603 


36 


64 


96 


13.6150(25) 


13.6141 


13.6476 


13.6155 


40 


64 


96 


14.8520(27) 


14.8492 


14.8757 


14.8499 


44 


64 


96 


16.0688(27) 


16.0699 


16.0918 


16.0701 


48 


64 


96 


17.2816(29) 


17.2804 


17.2993 


17.2804 



Table 9: Same as in table [71 but at (3 = 0.223952. 





L2 


Ls 




(/3), full 


F*(/3), 1-1. 


(/?), 2-1. 


24 


64 


96 


18.4131(26) 


18.4121 


18.4555 


18.4152 


28 


64 


96 


21.2414(27) 


21.2450 


21.2724 


21.2463 


32 


64 


96 


24.0310(27) 


24.0306 


24.0497 


24.0312 


36 


64 


96 


26.7859(28) 


26.7873 


26.8017 


26.7875 


40 


64 


96 


29.5271(30) 


29.5250 


29.5365 


29.5251 


44 


64 


96 


32.2449(31) 


32.2498 


32.2594 


32.2498 


48 


64 


96 


34.9623(33) 


34.9653 


34.9736 


34.9653 



Table 10: Same as in table [71 but at /3 = 0.226102. 
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Similar to the case of square interfaces, the numerical values of the full 
NG result and the two-loop approximation are very close down to \foL\ ~ 
1.8, i.e. they can not be discriminated, given the accuracy of our Monte 
Carlo results for the Ising model. 

In the range ^/oLx ^ 1.8, the Monte Carlo data perfectly agree with the 
full NG result and the two-loop approximation. Up to our largest values of 
\faL\^ there is a clear discrimination against the 1-loop approximation. 

For y/aLi < 1.8 the Monte Carlo data are more compatible with the 
2-loop approximation than with the full NG result. 

In the case of the Polyakov-loop correlator, for large distances of the 
Polyakov-loops, similar observations were made with respect to a finite tem- 
perature; i.e. for the direction with periodic boundary conditions [3-12]. 
Even for \/aLx close to ^/a/T^., the numerical data follow closely the 2-loop 
prediction. 

In our opinion, this behaviour at very small distances is a mere coinci- 
dence, related to the fact that the 2-loop approximation gives (by chance) 
the same critical exponent as the 2D Ising universality class and also the 
value for T^/^fa is very close to the one of the Z2 gauge theory. 



6 Universal amplitude ratios 

In this section we compute universal amplitude ratios of the interface tension 
and the correlation length. Generically, these amplitude ratios have the 
form: 

R = \\^^a{mtf , (45) 

where t is the reduced temperature. In particular, we shall consider the 
second moment correlation length in the high temperature and the low tem- 
perature phase and the exponential correlation length (i.e. the inverse of 
the lightest mass) in the low temperature phase. 

6.1 The second moment correlation length in the high tem- 
perature phase 

The second moment correlation length is defined by: 



= 4 , (46) 



where: 

/x = ^x2(soSx) (47) 
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Figure 4: In the two figures we give — aLiL2 as a function of \foL\ 

for /? = 0.223102 and fi = 0.223452. In all cases Lq = 96 and L2 = 64. Note 

that in the case of the Monte Carlo results the statistical error is smaller 

than the symbol (circle). The 1-loop, 2-loop and full NG predictions are 

given as solid black, red and blue lines, respectively. 
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Figure 5: Same as in the previous figure, but for /? = 0.223952 and (3 = 

0.226102. For /? = 0.226102, the 2-loop and full NG predictions fall (within 

the resolution of our plot) on top of each other. 
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and X is the magnetic susceptibility: 



X = ^(soSx) ■ (48) 

X 

The coefficients of the high temperature series for and x up to O{0^^) are 

reported in ref. [76]. Using these, we computed the coefficients for ^2nd- Quite 
a straightforward method is the so-called matching method (see e.g. [77]): 
We start from a power-law ansatz like: 

eL = /L,+(/3c -/?)-'" (49) 

[1 + a(/3, - I3f + 6(4 -0) + c{(5^ - f3f' + d(/3, - (3)'^ + ...], 

where we fix /?, = 0.22165455, u = 0.6302, 6 = 0.5174 and 7 = 1.2372, 
following the Monte Carlo results of ref. [68] . Now the ansatz (149p is Taylor- 
expanded in /?. The remaining free parameters of the ansatz (j49p can then 
be determined by matching the coefficients with those of the exact result for 
the HT-series for 6.2nd- our final result for the amplitude, we obtain: 

fL,+ = 0.038369(13) (50) 
X [1 + 1024 - 0.22165455) - 17 (1/ - 0.6302) - 0.0206 {9 - 0.5174)] . 

Combining the results of eqs. ()41|50p . we get: 

R+ = 0.3869(14) (51) 
+1250(/?, - 0.22165455) + 0.11{iy - 0.6302) - 0.017(6* - 0.5174) . 

Inserting the errors of the input parameters, we arrive at: 

R+ = 0.387(2) (52) 

as our final estimate. The most precise theoretical estimate given in the 
literature ii+ = 0.377(11) [74, 78] is fully consistent with ours. Similar to 
our analysis, this result is derived from HT-series of the correlation length 
and Monte Carlo data of the interface tension. 

Experimental measurements of this quantity are consistent with but less 
precise than our estimate; for instance, i?+ = 0.41(4) was obtained from the 
study of a cyclohexane-aniline mixture in [79] . 
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a 


^cxp 








0.23910 


0.055415 


1.296(3) 


0.0931(4) 


1.2335(15) 


0.0843(2) 


0.23142 


0.027601 


1.868(3)* 


0.0964(3) 


1.8045(21) 


0.0899(2) 


0.22750 


0.014740 


2.593(3)* 


0.0991(2) 


2.5114(31) 


0.0930(2) 


0.22600 


0.010228 


3.135(9) 


0.1005(6) 


3.0340(32) 


0.0942(2) 


0.22400 


0.004761 


4.64(3) 


0.1025(13) 


4.509(6) 


0.0968(3) 


0.22311 


0.002626 






6.093(9) 


0.0975(3) 








0.1084(11) 




0.1024(5) 



Table 11: The values for the interface tension are taken from our global fit, 
the values for the exponential correlation length are taken from ref. [81] and 
in the cases marked by * the results of ref. [81] and ref. [82] are averaged. 
The numbers for the second moment correlation length are all taken from 
ref. [80]. In the last row we give the extrapolation to the scaling limit. For 
details see the text. 

6.2 The correlation length in the low temperature phase 

Here we use the Monte Carlo results for the second moment correlation 
length ,^2nd obtained in ref. [80] and the exponential correlation length from 
refs. [81,82] to compute the combination a{t)^{t)'^ at finite values oft. With 
our present results for the interface tension, the error of the combination 
is completely dominated by the error of the correlation length. We fit the 
combination with the ansatz: 

ai^ = R^ + c , (53) 

using uj = 0.821(5) [68]. In the case of the second moment correlation length 
we take our final estimate from a fit that includes (3 < 0.2275. The depen- 
dence of the result on the value of uj is rather small and can be neglected 
compared with the statistical error of R-. In the case of the exponential 
correlation length our final estimate is obtained from data with /3 < 0.23142 
where, again, the error due to the uncertainty of uj can be ignored. The 
values of the interface tension, the correlation length and the amplitude 
combination i?_ as well as the final results of our fits are summarized in 
table dH 

For comparison, table [12] reports previous estimates obtained for the 
Ising model in refs. [56,66,74,81,83]. 

Also field-theoretical predictions for the universal constant R are given 
in the literature: 
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year 


authors(s) 


Ref. 


R- 


1992 


Klessinger and Miinster 


[83] 


0.090(3) 


1993 


Hasenbusch and Pinn 


[66] 


0.090(5) 


1996 


Zinn and Fisher 


[74] 


0.096(2)* 


1996 


Agostini et al. 


[81] 


0.1056(19) 


1997 


Hasenbusch and Pinn 


[56] 


0.1040(8)* 



Table 12: Comparison of a number of estimates for i?_ taken from the 
literature. The estimate of Zinn and Fisher is based on data of [66]. Agostini 
et al. used the true instead of the second moment correlation length. The 
result marked with a star refer to the second moment correlation length, 
while the others refer to the exponential one. 



From the e-expansion: Brezin and Feng 
e-expansion. Their result reads: 



computed R from the 



1 



2-K 



AirR 



/ 47 1 1 

1 - e h - ln(47r) 7 

I 54 2 ^ ^ 2 ' 
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+ 0(e3) , (54) 



with 7 = 0.5772.... Their numerical evaluation for e = 1 gives results 
in the range from R ~ 0.051 up to ~ 0.057. Note that this result 
deviates by about a factor of 1/2 from ours. 

From perturbation theory in 3D fixed: The study of interfaces using 
perturbation theory in three dimensions was pioneered by Miinster [85]. 
The starting point of this calculation is the classical solution (i.e. the 
configuration with minimal action) of a system with fixed boundary 
conditions in 3-direction. At one boundary, the field is fixed to the 
negative minimum, and at the other, to the positive minimum of the 
potential. Then, fluctuations around this classical solution are stud- 
ied. In the more recent paper [86], this analysis was extended to two 
loops. Their resullEI is: 



R 



''R 



with: 



37 
32 



-0.2002602 



(55) 



(56) 



'Note that they use the second moment definition of the correlation length. 
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and: 

(72; = -0.0076(8) . (57) 

Using Pade and Pade-Borel analysis, with = 14.3(1), they obtain: 

R = 0.1065(9) (58) 

as their final result. Note that the quoted error is dominated by the 
error of ujj. At this point, one should also note the principle problems 
related with the definition of 

ctS discussed in ref. [82]. 

Our value cr^^^p = 0.1084(11) corresponds to mo++/\/o' = 3.037(15). 
This value for the Z2 gauge theory might be compared with 4.718(43), 
4.329(41), 4.236(50) and 4.184(55) for the SU{2), SU{2>), SU{4) and SU{5) 
gauge theories in 2+1 dimensions, respectively [87]. Note that there are 
clear differences between the results for the different gauge groups, indicat- 
ing that the 0++ glueball state probes short distances, which can not be 
described by an effective string. This is in contrast to higher exited glueball 
states, where much less dependence on the gauge group is found [81,87]. 

In a similar way, we have analysed cr/T^, using the results of ref. [67] 
for T,. We arrive at a/T^ = 0.656(2) or TJ^/B^ = 1.235(2), where we have 
included data with (l/T^.) > 8 into the fit. Again the error is completely 
dominated by the error of T^. Our new result can be compared with our 
previous estimate T^j ^ = 1.2216(24) [67], (where we did not take into 
account scaling corrections, and had the interface tension only available up 
to l/T, = 12) and T./y^ = 1.17(10) [88]. 

Our result is clearly different from the value i/3/7r = 0.977... obtained 
from the effective string picture [54]. 

Finally we might compare with SU{N) gauge theories in 2+1 dimensions: 
In the literature one finds TJ^ = 1.12(1), 0.98(2) for AT = 2,3 [89] and 
TJ^ = 0.892(3), 0.879(3), 0.877(3) for A^ = 4,5,6 [90], respectively. The 
results for A" = 4, 5, 6 are read off from figure 4 of ref. [90]. Note that they 
are taken for Aj = 3 and no continuum extrapolation is performed. 

Similar to the case of mo++ / we observe a clear dependence on the 
gauge group. Nevertheless, it is quite remarkable that the simple string 
picture gives the correct value with less than 30% deviation. 

7 Conclusions 

In this paper, we have presented the results of an accurate numerical study 
of the interface free energy in the three-dimensional Ising model. The moti- 
vation for this study was twofold: 
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• To investigate the dynamics of a fluctuating interface with periodic 
boundary conditions only, avoiding non-trivial effects of Dirichlet bound- 
ary conditions. 

• To obtain high precision estimates of the interface/string tension in 
a large range of the inverse temperature, allowing us to compute the 
scaling limit of several universal amplitude ratios with high precision. 

To this end, we have determined the interface free energy by numerical 
integration of the interface energy over the inverse temperature (3. The 
interface energy is given by the difference of the internal energies in a system 
with periodic and a system with antiperiodic boundary conditions in one 
direction. The interface free energy at the starting point of the integration 
was measured using the boundary-flip algorithm. For our simulations we 
have used efficient combinations of cluster and multispin coded Metropolis 
updates. This approach has allowed us to strongly improve the precision 
of the numerical results, as compared to analogous studies presented in the 
literature. In particular, the large range of /^-values that we have studied 
gives us a good control over corrections to scaling (or, in the language of 
lattice gauge theory: finite a effects). It turns out that, in the regime 
not too close to the finite temperature transition, the interface free energy 
Fs{y/aLi, ^/aL2) approaches its continuum limit quite fast, characterized 
by a correction exponent uj' ~ 2. This observation might be explained by 
the fact that the effective interface model only assumes the restoration of 
the symmetries of the continuous space-time, which indeed comes with an 
exponent uj' ^ 2. 

The level of precision in the Monte Carlo data, as well as the accu- 
rate control of the systematic errors, enables us to clearly resolve the fine 
string-dynamics effects that we were seeking after: In a setting in which 
the possible distortions due to boundary effects are completely absent, and 
all systematic effects are under control, the Nambu-Goto model yields an 
accurate description of the data up to the second loop order (only). 

Comparing the 2-loop approximation with the full NG prediction one 
has to notice that down to rather small scales, such as ^/a^L ~ 1.8, the two 
string results can not be discriminated at the level of the accuracy of our 
Monte Carlo data. Nevertheless, one might interpret our Monte Carlo data 
as a confirmation of the full NG prediction in the sense that e.g. l/{aA)'^ 
corrections have indeed a small amplitude. 

Going close to the finite temperature transition, which occurs at y/aNt = 
0.810(2) as discussed in section [6l it does not come as too big a surprise that 
the data for the Ising model are not well fitted by the full NG prediction: 
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The full NG predicts mean-field behaviour of the transition, while following 
the Svctitsky-Yaffe conjecture [55] (numerically confirmed e.g. in [67]) it 
should be the behaviour of the 2D Ising universality class. 

This behaviour is quite similar to that of interfaces with a cylinder-like 
geometry, with periodic boundary conditions in the short direction (see e.g. 
ref. [12]), thus confirming the assumption that a common effective string 
description underlies the confining flux tube dynamics in different physical 
settings. 

These observations indicate that the all-order prediction of the Nambu- 
Goto effective string action (at least as it is treated in the approximation in 
which the role of the Liouville field is neglected) does not show a quantitative 
agreement with the data for an interface in the three-dimensional Ising model 
that goes beyond the two-loop approximation. This fact is — as discussed 
above — consistent with the Polchinski-Strominger model. 

The apparent failure of the two-loop prediction with respect to Dirichlet 
boundary conditions (see e.g. ref. [7]), as they arise via duality from the 
Polyakov-loop correlator requires further investigation. At least the present 
work confirms that the problem is intrinsically related with the boundary 
conditions. A possible explanation is that the Dirichlet boundary conditions 
probe short distance properties of the theory, which are not captured by the 
effective string model. 

Using our new data for the interface tension along with the analysis 
of the high temperature series expansion of the second moment correlation 
length, we have computed the universal amplitude ratio = 0.387(2). 
This estimate is more precise than any other theoretical estimate given in 
the literature. This estimate can be compared with experimental results e.g. 
for binary mixtures. In ref. [79] i?+ = 0.41(4) was found, which is consistent 
but less precise than our result. 

Finally, we have also updated the results for the mass mo++ / -v/o" of the 
0+"^ glueball in units of the square root of the string tension and the critical 
temperature of the deconfinement transition T^j ^fa. 

The comparison with other gauge theories in 2 -|- 1 dimensions shows a 
variation of these quantities by roughly 50%, indicating that these quanti- 
ties can not be described only by an effective string picture, but that also 
microscopic features of the gauge theory have to be taken into account. 
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A Integration schemes 

One of the simplest methods given in textbooks is the trapezoid rule: 

+ 0(iV-2), (A.l) 



f{x)dx = h 



+ h+h + ... + fN-l + \fN 



where h = [xjsj — xo)/N, fi = f{xi) and Xi = xq + ih. There are rules with 
faster convergence as ^ oo, like the well-known Simpson rule: 



Jxo 



f{x)dx = h 



1 4 2 4 2 



4 2 4 1 ■ 

+-^fN-3 + 3/^-2 + 3/^-1 + 3/^ 



+ 0(iV-^) . (A.2) 



However, the disadvantage of the Simpson rule is that fi is not constant in 
the middle of the interval: this implies a loss of precision in the final result. 
Furthermore, N has to be even. 

A better-suited rule that avoids these problems is given by (see e.g. 
eq. (4.1.14) in ref. [91]): 



r 

Jxo 



f{x)dx = h 



3 7 23 

g/0+g/l + 24/2 + /3 + /4 + ... 

23 7 3 



+ 0{N-^). (A.3) 



A similar and maybe slightly better rule (see e.g. eq. (35) in ref. [92]) is 
given by: 



/ f{x)d 



x = h 



17 59 43 ,49 

4^^o + 48^i + 48^2 + ^/3 + /4 + 



,f I 49 43 59 17 



+ 0(Ar-4).(A.4) 
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